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1. INTRODUCTION 


Experiments have been performed to investigate the convective heat 
transfer for transverse flow over solid cylinders located concentrically 
within perforated cylindrical shrouds. This work was motivated by the need 
for convective heat transfer data applicable to these configurations. Such 
data could be used in the development of response models for platinum resis- 
tance thermometers (prt's) which have similar shrouded cylinderical configu- 
rations. One of the important parameters in the response modeling of prt's 
is the convective heat transfer coefficient to the sensor (cylinder) [1].* 

Utilization of perforated shrouds on bare cylinders can serve two 
purposes. First, they act as radiation shields. Perforated shrouds have 
also been shown to be effective in reducing vibration of cylindrical geom- 
etries caused by vortex shedding [2]. In an effort to develop ways of 
suppressing wind-induced oscillations of tall smokestacks. Price [3] has 
concluded that the enclosure of a cylinder within a concentric perforated 
shroud was an effective method of vibration reduction. He has also found 
that the drag coefficient of shrouded cylinders was approximately 0.6, and 
it was relatively immune to variations in Reynolds number. Walshe [4] has 
also reported that perforated shrouds could stabilize flow-induced oscil- 
lations around cylinders. While reducing vibrations of cylinders, shrouds 
tend to complicate the convective heat transfer to cylinders. The influence 
of these perforated shrouds on the convective heat transfer to the cylinder 
located concentrically inside them is the subject of this study. 

No definitive heat transfer correlations for transverse flow over cyl- 
inders located concentrically inside perforated shrouds exist. Warshawsky 
[5] has addressed this issue, and has recommended using the shroud temper- 


*Number in brackets represents references. 



ature as the free stream fluid temperature and estimating a fluid velocity 
based on perforation size and internal flow impedance, in order to calculate 

the heat transfer coefficient for the inner cylinder. He then recommended 

using conventional bare cylinder heat transfer correlations for this ideal- 
ized case. Those approximations have not been substantiated, and that was 
an additional consideration in this experimental investigation. 

The convective flow through the shroud represents a formidable theo- 
retical problem. Since patterns of circular holes are used typically to 
allow fluid flow to occur in the annular space between the shroud and the 
inner cylinder, the flow is periodic in the spanwise direction and three- 
dimensional. Furthermore, the presence of the ventilating holes provides 
pressure relief around the outer cylinder so that classical solid cylinder 
pressure distributions around the shroud do not occur. In fact, it is the 

balancing of the two pressure distributions — one due to classical flow 

around cylinders and the other due to suction and blowing through the venti- 
lation region -- that controls the two coupled flow fields. These flow 
fields are influenced by the following parameters: 

1. Reynolds ntnrber 

2. Degree of shroud ventilation 

3. Shroud orientation with respect to the flow 

4. The size of the flow passage between the shroud and the inner 
cyl inder 

5. Free-stream turbulence intensity 

6. Aspect ratio — length to diameter ratio of the sensor assembly 

7. Perforation hole size and shape 

8. Perforation hole pattern 
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9. Thickness of shroud 
10. Surface finishes 

We can define a shroud ventilation factor, t|», as the ratio of the open area 
to the total surface area of the perforated portion of the shroud. The flow 
between the shroud and cylinder can be characterized by a radius ratio, <j>, 
which is the ratio of shroud's inside radius to the radius of the inner 
solid cylinder, and serves as a nondimensional parameter representing the 
clearance between the two tubes. 

It was not possible to investigate all the parameters defined here. 
However, the influence of Reynolds number, ventilation factor, and shroud 
orientation with respect to flow were studied. Some of the results were 
presented in an earlier publication [6]. In the continuation of this inves- 
tigation it was decided to include also the influence of radius ratio, be- 
cause it appeared to be a very important parameter. The present work pro- 
vides comprehensive results obtained throughout this study. 

In order to place the shrouded cylinder work in perspective, a brief 
discussion of earlier bare cylinder experiments is required. Analytical 
solutions of the convective heat transfer coefficient for different solid 
geometries require the solution of governing conservation equations (mass, 
momentum, and energy), and are restricted to relatively simple geometries. 
Convective heat transfer coefficients for most geometries, including cylin- 
ders, are evaluated from empirical relations determined by a combination of 
dimensional analysis and experiments. 

The three most important nondimensional i zed parameters for convective 
heat transfer correlations are: Nusselt number, Nu, which is the basic 

dimensionless convective heat transfer coefficient (Nu = lu/k), Reynolds 
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number. Re, which is the ratio of inertial to viscous forces of the fluid 

Ujj, 

(Re = ) , and Prandtl nunber, Pr, which is the ratio of mechanical to 

V 

uC 

thermal diffusion (Pr = l = ^.). 

k a 

Numerous correlations for forced convection heat transfer for cylin- 
ders in cross flow exist. A summary of recent work on the mean heat trans- 
fer rates from circular cylinders is contained in a paper by Morgan [7]. 
Morgan has reported that the Nusselt number for a given Reynolds number 
varies by up to 29% for various experimental data and as much as 46% for 
various correlations. These large variations are caused by such factors as: 
experimental uncertainty, and nonsimilarity of experimental conditions 
(different surface roughness, free-stream turbulence intensities, etc.). 

Measurements of forced convective heat transfer coefficients for cyl- 
inders can be found in papers by Hi 1 pert [8] and Schmidt and Wenner [9]. 
Hilpert performed measurements on circular cylinders in a cross flow of air 
covering a wide range of Reynolds numbers. Schmidt and Wenner [9] have 
measured the variation of local heat transfer coefficients over the surface 
of cylinders. Giedt [10] has also measured the variation of local heat 
transfer coefficients over non-i sothermal cylinders. 

The correlation used in this investigation for the mean convective 
heat transfer rates from cylinders is that proposed by Morgan [7] which is 
the corrected form of Hi 1 pert's correlation [8]. 

2. EXPERIMENTAL SET-UP 

For the purpose of clarity, a typical industrial platinum resistance 

thermometer is shown schematically in Figure 1. It can be seen that the tip 
of the sensor consists of two concentric cylinders — an inner cylinder 
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Fig. 1. Schematic diagram of a typical industrial platinum resistance thermometer. 







encapsulating the resistance element, and an outer perforated cylindrical 
shroud fitted over the inner cylinder. That basic geometry was used in 
designing the present experiments. 

All of the experiments in this investigation were conducted in the 
NASA, Langley Research Center Instrument Research Division's low speed, open 
circuit wind tunnel. The tunnel was manufactured by Astro Dynamics and has 
a test section which is 30 cm x 43 cm. The facility could be operated with 
flow speeds between 3.5 and 80 m/s. The wind tunnel motor could be operated 
at both low and high speed ranges. When the wind tunnel motor was set at 
the low speed range, flow velocities between 3.5 and 40 m/s could be 
obtained by controlling vane openings. All of the experiments in this 
investigation were conducted with the motor at low speed. The variation of 
the free-stream turbulence intensity with free-stream velocity is shown in 
Figure 2 for the motor set at the low speed range. The data were obtained 
using hot wire anemometry. The high turbulence intensities at the lowest 
speeds may be related to excessive instrumentation noise. 

The experimental apparatus consisted of an inner cylinder, referred to 
as the "test cylinder" or "bare cylinder," and a set of perforated shrouds 
with different ventilation ratios and different diameters. A photograph of 
the test cylinder and some of the shrouds used in the experiment is shown in 
Figure 3. This figure also includes the end mounts and struts used for 
installing the experimental apparatus in the wind tunnel. A photograph of a 
shroud fitted over the inner cylinder is shown in Figure 4. A photograph of 
the experimental apparatus installed in the wind tunnel is shown in Figure 
5. The wind tunnel and the supporting instrumentation used in this investi- 
gation are shown in Figure 6. 
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Fig. 2. Free stream turbulence intensity of wind tunnel facility at low speeds, courtesy of 
Debra L. Carraway, NASA Langley Research Center. 
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Fig. 4. A perforated shroud fitted over the bare cylinder. 
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Fig. 5. Experimental apparatus installed in the wind tunnel 
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and related instrumentation. 


2.1 Test Cylinder 


A brass tube with an outside diameter of 7.9 mm and a wail thickness of 
1.1 mm was used as the test cylinder throughout the experiments. Since the 
local heat transfer coefficient on a circular cylinder varies considerably 
across the surface of the cylinder [9], a good thermal conductor is required 
to approximate an isothermal surface. Therefore,, brass, a relatively high 
thermal conductivity alloy, was used in order to minimize surface temper- 
ature variation (straight copper tubes with the desired outside diameter 
were not available). However, local temperature variations were assumed 
throughout the experiment. 

The test cylinder had an aspect ratio (length to diameter ratio) of 26. 
It was decided to use the central 76.2 mm of the length of the cylinder, 
where longitudinal variation of temperature is negligible (refer to Appendix 
A), for the experimental determination of the convective heat transfer co- 
efficient. A sectional view of the test cylinder with the specified 76.2 mm 
instrumented section is shown in Figure 7. 

Five 0.25 mm diameter kapton insulated chromel/alumel thermocouples 
were mounted near the surface of the test cylinder for surface temperature 
measurement. The thermocouples were staggered 12.7 mm on either side of tne 
test cylinder symmetry plane, and were located at angles of 0° (forward 
stagnation line), 15°, 75°, 120°, and 300° as shown in Figure 8. The therm- 
ocouple wires were drawn through the central region of the test cylinder, 
then out to the surface through 0.64 mm diameter holes, and soldered in 
place. The thermocouple junctions could not be installed exactly at the 
surface. It was estimated that the thermocouple junctions were located to 
within one wire diameter from the tube surface. Calculations showed that a 
position error equal to one wire diameter into the tube, would produce a 
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measured surface temperature error of less than 0.3°C at the steady state 
heating rate used throughout the experiment (1040 w/m 2 ). That estimate was 
conservative since it was assumed that the thermal environment was con- 
trolled by a low thermal conductivity potting compound, rather than by the 
brass. The potting compound was used to fill the tube and will be described 
later. The general installation of a thermocouple in the test cylinder is 
also shown in Figure 7. During preliminary testings it was found that the 
thermocouple junction located at 15° was not close to the surface; there- 
fore, measurements at that location were ignored. 

An 18 gauge (1 mm diameter) chrome! wire located on the axis of the 
test cylinder served as the heating element. The centrality of the heater 
wire was very important in producing uniform heating across the test cyl- 
inders. Two 0.08 mm diameter teflon insulated chromel wires were spot 
welded to the heater wire for measurement of the voltage drop across the 
central 76.2 mm of the test cylinder assembly as shown in Figure 7. 
Simultaneous measurements of the voltage drop across the central portion of 
the heater wire, and the current supplied to the heater wire were used to 
determine the steady state energy flow per unit length in the central 
portion of the test cylinder. 

The heating wire, voltage leads, and thermocouples were potted in place 
using EpoxyLite #205 (manufactured by EpoxyLite Corp.)* which had a thermal 
conductivity of 0.52 W/mK for the temperature range of this experiment. The 
following procedure was performed for applying the potting compound: 


*Certain commercial equipment, instruments, or materials are identified in 
this paper in order to adequately specify the experimental procedure. Such 
identifications do not imply recommendation, nor does it imply that the 
material or equipment identified are necessarily the best available for the 
purpose. 
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1. The test cylinder was set vertically. 

2. The heating wire was set at the central axis of the tube, held in 
position by applying tension. 

3. The top and bottom part of the tube were capped. 

4. The potting compound was pimped into the tube from the bottom part 
using a syringe. 

5. The potting compound was cured at room temperature for approxi- 
mately 12 hours. 

After the completed test cylinder was cured, the assembly was x-rayed to 
verify the position of the heater wire. It was found that the heater wire 
was centered to within the resolution of the x-ray picture. 

2.2 Perforated Shrouds 

Perforated shrouds employed on industrial platinum resistance thermom- 
eters typically have 4 to 6 rows of staggered circular holes, with venti- 
lation factors of .15 to .20, and radius ratios between 1 and 2. Three sets 
of stainless steel shrouds were chosen for this investigation that would 
result in radius ratios of 1.1, 1.4, and 2.1 when fitted over the test cyl- 
inder. Stainless steel tubes were selected because of ease of punching the 


required holes 
sions of each 

on their surfaces 
set of shrouds are 

without deforming the tubes, 
listed in Table 1. Each set 

The dimen- 
consisted of 

Table 1. Dimensions of test shrouds. 
Nominal 

Actual 

Radius Ratio 

O.D. 

Wall Thickness 

Radius Ratio 

1.1 

9.5 mm (3/8") 

.25 mm (.010") 

1.14 

1.4 

12.7 mm (1/2") 

.89 mm (.035") 

1.38 

2.1 

19.05 mm (3/4") 

1.24 mm (.049") 

2.09 
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three perforated shrouds with nominal ventilation factors of 0.09, 0.18, and 
0.27. The set of perforated shrouds with radius ratio of 1.1 also included 
a tube with a ventilation factor of 0.36. Six rows of circular holes were 
drilled in the stainless steel tubes to produce the various ventilation 
factors. The holes were staggered in successive rows which were located at 
60° intervals around the tube circumference. The ventilation factors were 
controlled by using different hole diameters and hole spacings. Details 
about the hole size, hole spacing, and total nunber of holes used on each 


perforated shroud 

are provided in 

Table 

2. 




Table 2. 

Oesign data for perforated shrouds. 




Tube No. 

Radius 

Ratio 

* 

Nominal 
Ventil ation 
Factor 

* 

Hole 

Diameter 

(mm) 

Hole 

Spacing 

(mm) 

Total 
Number 
of Holes 

Actual 

Ventilation 

Factor 

1 

1.1 

.09 

3.2 

(1/8") 

18.5 

45 

.090 

2 

1.1 

.18 

3.2 


9.0 

87 

.180 

3 

1.1 

.27 

3.3 

(.130") 

6.4 

117 

.263 

4 

1.1 

.36 

4.0 

(5/32") 

6.9 

111 

.359 

5 

1.4 

.09 

4.8 

(3/16") 

30.6 

27 

.091 

6 

1.4 

.18 

4.8 


15.3 

51 

.173 

7 

1.4 

.27 

4.8 


10.2 

75 

.254 

8 

2.1 

.09 

6.3 

(1/4") 

40.2 

21 

.083 

9 

2.1 

.18 

6.3 


17.2 

45 

.178 

10 

2.1 

.27 

6.3 


12.1 

63 

.250 


The shroud temperature was required for calculation of radiation heat 
flux correction. Two 0.25 mm diameter chromel/alumel thermocouples were 
mounted on each shroud to measure characteristic shroud temperatures. The 
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thermocouples were attached to the outer surface of each shroud at nominal 
distances of 6.0 cm from one end of the shroud, separated by an angle of 
approximately 180°. That location was considered far enough from the end 
mount to minimize end effects, and also outside of the central instrumented 
section to minimize flow distortions. The thermal junctions of these therm- 
ocouples were attached to the shrouds by spot-welding the individual wires, 
approximately 1.6 rrm apart, to the outer surface of the shrouds. This 
installation, which causes the metal of the shroud to become involved in the 
thermoelectric circuit, has proven to be satisfactory when no temperature 
gradients exist within the (1.6 mm) region separating the two wires [11]. 
Only one thermocouple was attached to each of the perforated shrouds with a 
radius ratio of 1.1. 


2.3 Related Instrumentation 

The supporting instrumentation used in this experimental investigation 
is described in this section. A Hewlett-Packard DC power supply (model 
6267B, current range 0-10 amps) was used to provide power to the heater 
wire. For tests involving the perforated shrouds with radius ratios of 1.4 
and 2.1, two Trygon Electronics DC power supplies (model HR60-5B, current 
range 0-5 amps) were used in parallel, to supply the required 5.5 amps 
current to the heater wire. Voltage and current applied to the heater were 
measured with a Dymec integrating digital voltmeter (model 2401 C) , and a 
Fluke digital multimeter (model 8120A), respectively. Thermocouple temper- 
ature measurements were made using Joseph Kaye and Company, Inc. ice point 
reference standards, and the Dymec voltmeter. Wind tunnel dynamic pressure 
measurements were made using a Datametrics Barocel differential pressure 
transducer (model 570) and electronic manometer (model 1174). 
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3. EXPERIMENTAL PROCEDURE AND ANALYSIS 
The bare cylinder correlations of Hi 1 pert [8] have been studied by 
Morgan [7] and corrected using improved thermophysical property data. The 
correlation is in the form: 


Nu m = D 2 Re n i (1) 

where Q? and ni are the correlation constants, and are contained in Table 3 
along with Hilpert's original correlation constants. 


Table 3. 

Comparative correlation constants for 
to cylinders (for use with Equation 1) 

convective 

heat transfer 


Hilpert's Constants 



Morgan' 

's Constants 


Re 



Re 



From 

To D 2 

n i 

From 

To 

0 2 

"l 

1 

4 .891 

.330 

1 

4 

* 

* 

4 

40 .821 

.385 

4 

35 

.795 

.384 

40 

4X10 3 .615 

.466 

35 

5xl0 3 

.583 

.471 

4x103 

4x10** .174 

.618 

5xl0 3 

5xl0 4 

.148 

.633 

4X10 4 

4xl0 5 .0239 

.805 

5x1 0 4 

2.3xl0 5 

.0208 

.814 

*Too 

few data. 







Both Nusselt and Reynolds numbers are based on the cylinder diameter. 
It should be noted that all the correlations reported here are for the 
special case of air, but could be extended to other fluids by multiplying 
the Nusselt number by the ratio of the Prandtl number of the fluid to the 




Prandtl number of air raised to an appropriate power (usually 1/3). 

3.1 Test Procedure 

It was decided to perform the experiments for two sets of Reynolds 
nunbers given in Table 3. The Reynolds nunber ranges of 35 to 5,000 and 
5,000 to 50,000 were selected. The nominal Reynolds numbers used throughout 
the experiments were 1800, 2700, 4000, 6000, 9000, and 14000. This Reynolds 
number range was representative of Reynolds numbers that would occur on 
shrouded prt's located in the stagnation region of the National Transonic 
Facility at NASA Langley Research Center. Since this investigation was 
motivated by a need for improved response modeling of platinum resistance 
thermometers for that facility, that Reynolds number range was of primary 
interest. 

It was desirable to maintain a minimum temperature difference of 5°C 
between the surface of the test cylinder and free stream. A simple one 
dimensional -- heat transfer analysis (neglecting axial variation) for the 
multilayered cylinder was performed, and it was found that a heater current 
of 5.5 amps would be sufficient to maintain the desired 5°C temperature 
differential at the highest Reynolds number used in this experiment. This 
current would correspond to a 1040 w/m 2 heat flux based on the surface area 
of the instrumented section of the test cylinder. This heat flux value was 
used throughout the experiments. 

The first step in the experimental investigation consisted of measuring 
convective heat transfer coefficients for the test cylinder in normal flow 
over the Reynolds nimber range 1000 to 20000. The results were compared 
with Morgan's correlations [7] to validate the experimental procedure. 

Then, the experiments were repeated with the perforated shrouds fitted over 
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the test cylinder. 

The variation of convective heat transfer with Reynolds nunber, venti- 
lation factor and radius ratio of shrouds was studied. Originally, shroud 
orientation with respect to flow was not considered to be a primary param- 
eter in this study, but it became apparent that this effect could not be 
ignored. Referring to Figure 9, it can be seen that a variety of flow 
fields can be created in the annular region by rotating the shroud about its 
axis. Defining the reference orientation as 0 = 0° (case a), it can be seen 
that the two coupled flow fields can be altered continuously by rotating the 
shroud through angles of up to 30° from the reference orientation. At that 
point syrranetry becomes a factor. In order to assess the shroud orientation 
variable, the three configurations tested are described briefly in Table 4. 
For shroud cases with radius ratios of 1.4 and 2.1, only configurations a 
and c ( 0 = 0°, 0 = 30°) were tested. 


Table 4. Shroud 

hole configurations 

(Refer to Figure 9) . 

Designation 

Flow Angle 

Comments 

a 

o 

O 

II 

CD 

Symmetric, maximun pressure relief 

b 

0 = 15° 

Asymmetric, intermediate flow 

c 

o 

O 

CO 

11 

CD 

Symmetric, maxinuin pressure difference 


3.2 Surface Temperature Variation 

Local heat transfer coefficient on a circular cylinder varies consider- 
ably across the surface of the cylinder [9]. This would in turn cause cir- 
cumferential variation of surface temperature around the cylinder. It was 
desirable to consider average heat transfer coefficients in this investi- 
gation. Therefore, an average surface temperature was required in order to 
properly estimate the average heat transfer coefficient. One major concern 
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in this analysis was whether the four surface temperature measurements 
around the cylinder were sufficient to yield an accurate average surface 
temperature. However, it was believed that the circumferential variation of 
surface temperature and heat transfer coefficient would be more pronounced 
for the bare cylinder case compared to the shrouded cylinder case. The 
shrouded cylinder case was believed to have a much stronger viscous cou- 
pling, a smeared stagnation region, and with flow separation occurring much 
further downstream compared to the bare cylinder case. 

Even though a relatively high thermal conductivity tube (brass) was 
utilized to prevent significant surface temperature variation, further 
analytical investigation was required to determine maximum surface tempera- 
ture variation for the limiting bare cylinder case. 

Influence of the circumferential variation of heat transfer coefficient 
on the surface temperature of the test cylinder was studied using a two 
dimensional finite volume nunerical model of the test cylinder. A listing 
of the program is presented in Appendix B. Schmidt and Wenner's [9] data 
for the circumferential variation of Nusselt nimber, for a Reynolds number 
of 39,800 (representative of the highest Reynolds number used in this exper- 
imental investigation), shown in Figure 10, was used in this study. It was 
found that the surface temperature variation would be negligible, with a 
maximum variation of 0.03°C under the simulated test conditions. Based on 
these results, it was concluded that using four surface temperature measure- 
ments around the cylinder, for calculating the average surface temperature 
was sufficient. 


3.3 Determination of Convective Heat Transfer Coefficient 
The electrical heat flux applied to the test cylinder can be related to 
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(1) Re = 39800 

(2) Re= 101300 

(3) Re = 170000 

(4) Re = 257000 



(5) Re = 426000 


Fig. 10. Circumferential variation of convective heat transfer around a cylinder, after Schmidt 
and Wenner [9] . 





the desired average heat transfer coefficient, Tf, through Newton's law of 
cool ing: 


h = (q"-q;;)/(T s -Tj (2) 

where: 

q" = applied heat flux 

qj^ = radiation heat flux correction 

T $ = average test cylinder surface temperature 

T ro = free- stream temperature 

The simultaneous measurement of the voltage drop across the heater (V), and 
the current applied to the heater (I), could be used to calculate the 
steady-state energy flow rate in the central, instrumented section of the 
test cylinder (Q = V I). The energy flow rate could be converted directly 
to an energy flux, q", through the outside surface of the central section of 
the test cylinder, which is given by: 

q" = 0 /tt D s L (3) 

where D g = diameter of test cylinder 

L = length of the central instrumented section of test cylinder, 

76.2 mm 

The average surface temperature of the test cylinder was determined using 
the surface mounted thermocouples. No attempts were made to correct the 
measured surface temperatures for temperature field distortions caused by 
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the presence of the thermocouples. Free-stream temperatures were measured 
with a 0.25 mm diameter kapton insulated chromel/al umel thermocouple located 
in the wind tunnel stream, 2.5 cm upstream and 7.6 cm below the axis of the 
test cylinder. Upper and lower wind tunnel surface temperatures were also 
measured using 0.25 mm diameter kapton insulated chromel/al umel thermo- 
couples in order to check for flow stratification. Maximum temperature 
differences of 1°C were observed between the upper and lower walls. The 
stratification effects were neglected throughout the analysis. However, the 
wind tunnel wall temperature measurements were used in the calculation of a 
radiation heat flux correction for the bare cylinder cases. 

For the case of the bare cylinder, gray body, non-participating medium 
assumptions were used to estimate the radiation exchange between the test 
cylinder and the wind tunnel walls. The corresponding radiation heat flux 
correction was: 


q R ' r^T [ J > ■ a T sl < 4) 

where: 

e = total emissivity of test cylinder 
ct = Stephan-Boltzman constant 
= radiosity of the test cylinder 

The derivation of the above equation from the basic governing laws of radi- 
ation heat transfer, and the calculation of the test cylinders radiosity are 
contained in Appendix C. The outside surface of the test cylinder was 
coated with silver acrylic lacquer (with a measured total emissivity of 
0.36) so that the radiation heat exchange with the surroundings could be 
estimated accurately. 
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Estimation of radiation correction for the case of shrouded cylinders 
proved to be more difficult, since the test cylinder would exchange radi- 
ation with the perforated shroud, and the wind tunnel walls through the 
perforations on the shroud. The shroud would also exchange radiation with 
the wind tunnel walls. Exact radiation modeling of this problem was formi- 
dable and excessively time consuming. Consequently, two different approxi- 
mate methods were used: 

1. Neglect the presence of the shroud, and consider radiation exchange 
only between the bare cylinder and the wind tunnel walls. 

2. Consider the shroud to be unperforated, thus assuming radiation 
exchange between two infinitely long concentric cylinders with gray 
surfaces through non-participating medium. 

It was found that both approximations produced comparable results. The 
maximum radiation heat flux (at the lowest Reynolds nimber) amounted to less 
than 5% of the applied I 2 R heat flux. It was decided to use the second 
approximate method for calculating the radiation correction. The corre- 
sponding correction was calculated [12] according to: 


- T sh 4 > 

1 + hi + °s (1 ' e sh ) 

E IO sh c sh 


where: 

e gh = total amissivity of shroud 

ID = shrouds inside diameter 
sh 

T s ^ = average shroud temperature 


(5) 
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The heat transfer coefficients obtained from Equation 2 were converted 


into Nusselt number data through: 


Nu 



( 6 ) 


where k is the thermal conductivity of airstream. All the thermophysical 
properties for the airstream were evaluated using dry air tables [13] at the 
mean temperature, T , as suggested by Eckert [14]: 


T + T 


T = 
m 


(7) 


The experimental Nusselt numbers were compared with Morgan's correlations 
[7] for cylinders in normal flow as given by Equation 1. The Reynolds 
numbers were calculated using 


Re = i (8) 

V 

where: 

U = free-stream velocity 
v = kinematic viscosity of airstream 

The free-stream velocities were calculated using Bernoulli’s equation: 


28 



(9) 


where: 


U = 



p = airstream density 
= dynamic pressure 

Dynamic pressures were measured using a pitot-static tube located in the 
* 

test section of the wind tunnel, 30 cm ahead of the test cylinder. The 
free-stream density was calculated using the ideal gas equation: 


where: 



( 10 ) 


R = gas constant for airstream 

P = static pressure 

Factors which complicate experimental heat transfer investigations are: 
wind tunnel blockage, free-stream turbulence intensity, and natural convec- 
tion. These factors have been examined, and are addressed subsequently. 

3.3.1 Flow Blockage Effects 

The flow interference by the model and its supports in wind tunnel 
testing causes the velocity at the location of the model to be higher than 
the undisturbed wind tunnel speed. Two of the most important components of 
this interference are: 

1. Solid blockage which causes an increase in the axial velocity at 
the model location due to a reduction in the cross sectional flow 
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area of the wind tunnel. Growth of the wind tunnel wall boundary 
layer also causes small increases in the free stream velocity. 

2. Wake blockage which is similar to solid blockage, but is caused by 
the reduction of speed in the wake of the model. That reduced 
speed wake causes the mass flow rate to increase outside of the 
wake to maintain conservation of mass. 

Blockage is usually expressed in terms of a blockage factor, X, defined by 
[15]: 


Up = U (1 + X) 


( 11 ) 


where U is the undisturbed tunnel speed and Up is the speed in the tunnel 
at the position of the model. Different forms of solid and wake blockage 
corrections exist [7], The solid blockage fraction, X ss , for a cylinder 
spanning the width of a closed rectangular tunnel, and located at a dis- 
tance, $, from the tunnel axis can be estimated [16] by: 


X = 0.822 (1 + — — ) (V (12) 

ss 4 e 2 e 

where e represents the height of the tunnel. The corresponding wake block- 
age factor, X w , can be estimated [17] by: 


X 


w 



(13) 


where is the drag coefficient. It is usual to assume that the correc- 
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tions are independent and simply additive [15], so the total blockage 
factor, X, for a cylinder with length, £, is: 

X = - {0.822 (1 + IZlf. ®L) (£) 1 2 + I (£) C n } (14) 

b 4 e 2 e 4 e u 

where 


b = width of tunnel 

The total blockage for the test cylinder and its cylindrical supports was 
calculated using the above equation. It was found that the maximum total 
blockage correction would amount to free stream velocity changes of approxi- 
mately one percent. That is, the local free stream velocity. Up, was: 

U F = 1.01 U (15) 

Even though this correction is very small, it was applied to all the free- 
stream velocity measurements. 

3.3.2 Free-Stream Turbulence Effects 

Some of the discrepancies observed among previous experimental results 
on heat transfer to cylinders can be attributed to the difference in free- 
stream turbulence intensities in the different investigations [7]. 

Increasing the intensity of free-stream turbulence will produce the follow- 
ing results [18]: 

1. Earlier transition to turbulence in the boundary layer, and hence 

an increase in heat transfer rate which is characteristic of a 
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turbulent flow as compared with a laminar boundary layer. 

2. Increase of heat transfer in the laminar portion of the boundary 
layer. At large turbulence intensities, a system of stationary, 
counter rotating vortices are set up in the boundary layer of the 
cylinder. Those vortices are the result of non-linear amplifi- 
cation of Taylor-Gortler instability modes in the transition 
process and produce a highly three-dimensional flow over the cyl- 
inder. The three-dimensional boundary layer flow increases the 
heat transfer rate. 

Kaylar [19] has investigated the effect of turbulence intensity on the 
transfer of heat from a circular cylinder. He has found that turbulence 
intensities below 1% have negligible effect on the heat transfer rate, as 
shown in Figure 11. The free-stream turbulence intensity of the wind tunnel 
in which this investigation was conducted was well below 1% as seen in 
Figure 2. It has been estimated that Hilpert's experimental data, on which 
Morgan's correlations are based, were obtained at free-stream turbulence 
intensity levels of 0.9% [20]. Therefore, it was concluded that the differ- 
ence in the free stream turbulence intensity of this experimental investi- 
gation and that of Hilpert's would not cause significant discrepancies in 
the heat transfer coefficients. 

3.3.3 Natural Convection Effects 

At high Reynolds numbers, heat transfer occurs mainly by the process of 
forced convection, but as the Reynolds number decreases, the contribution 
due to natural convection becomes significant. It was decided to investi- 
gate whether natural convection had a significant contribution to the heat 
transfer process at the lowest Reynolds number used in this experiment 
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Fig. 11. Effect of free stream turbulence intensity on the rate of heat transfer in the 
stagnation region of a cylinder, after L. Kaylar [19]. 


CJ 

OJ 



(1800). The nominal temperature difference between test cylinder's surface 
and free stream was about 15°C at this condition. This temperature gradient 
would yield a Nusselt number based on pure natural convection of about 3. 

The corresponding pure forced convection Nusselt number would be approxi- 
mately 20. 

Much confusion has existed over calculating the Nusselt number when 
both natural and forced modes of convection are present. Morgan [7] has 
recommended using an effective Reynolds number for calculating total heat 
transfer by Equation 1. This effective Reynolds number is given by: 

Re g ^ = yj( Re* + Re cos 5) 2 + (Re sin £) 2 (16) 

where: 

5 = angle between direction of flow and the vertical direction of 

natural convection flow 

Re* = equivalent Reynolds number for natural convection 
The equivalent Reynolds number is found by equating the Nusselt number for 
natural and forced flows. For natural flows, the Nusselt number is calcu- 
lated [7] by: 


Nu = B 2 (Gr.Pr) mi (17) 

where: 

Gr = Grashof number 

Pr = Prandtl number 

B 1 , m 2 = constants depending upon the magnitude of GrPr 
This yields: 
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Re* = 


l/r»i mj/nj 
Gr 


(18) 


R m i 

£L Pr 
D 2 


For the case of lowest flow speeds, Re = 1800, the equivalent Reynolds 
number is 26, having no significant effect on the effective Reynolds nunber. 
Thus, it was concluded that the effect of natural convection is negligible 
throughout the experiments. 


3.4 Error Analysis 

Measurement of any physical quantity involves some degree of error. 

The measured quantity is always disturbed by the act of measurement, which 
makes a perfect measurement practically impossible. Furthermore, all the 
measuring sensors used for typical experimental investigations have some 
degree of imprecision compared to the standard sensors — the most accurate 
sensors employed by national laboratories. One major concern in any experi- 
mental investigation is the estimation of the accuracy of the individual 
measurements. That is, the estimation of how close the measured values are 
to the "true" values. The imprecision of any measurement is estimated by 
using statistical analysis of the sensors calibration data, or manufac- 
turer's accuracy statements, or rough estimates based on the experience of 
the investigator. 

The overall accuracy of any physical quantity computed using a few 
independent measurements can be estimated using the root-sum-square method 
[21]. That is, if N is a function of n independent measurements, u 2 , u 2 . 



N = f (ui, u 2 , ..., u n ) 


(19) 
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where the independent measurements have imprecisions of ±AUj, ±au 2 . 


±Au n , respectively, then the overall accuracy of N is: 



where represents partial differentiation. 

3u . 

The overall accuracy of the Reynolds numbers and Nusselt numbers 
obtained in this experimental investigation were estimated using the root- 
sum-square method. In order to accomplish this all the dependent variables 
used in the calculation of Nusselt number and Reynolds number had to be 
expressed in terms of the independent measured quantities. The expression 
for the Nusselt number was (using Equations 2 and 6): 

Nu = LI (21a) 

TT L k (T - T ) 
s » 


where the effect of radiation exchange has been neglected for simplicity of 
analysis. The expression for Reynolds number was (using Equations 8, 9, 
10 ): 



The thermal conductivity and viscosity were not independent measurements. 
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and were calculated as a function of mean temperature, T . Therefore, they 
had to be expressed in terms of the independent measurements, T g and T^. 

Dry air data for these thermophysical properties [13] were curve fitted for 
the temperature range of 283 to 313 K: 


k = 3.92 x 10- 3 + 3.7 x 10- 5 

< T s 

8 

t— 

+ 

w/mK 

(22a) 

y = 4.39 x 10- 6 + 2.35 x 10- 8 

< T s 

+ 

— 1 
8_ ^ 

N.S/m 2 

(22b) 


The corresponding standard deviations of curve fitted data from the dry air 
tables were 1.68% and .01% for thermal conductivity and viscosity, respec- 
tively. These expressions were substituted in Equation 21a and 21b, and 
root-sum-square method was applied to the set of experimental data for the 
bare cylinders. The following imprecisions for the independent measurements 
were used: 

Temperatures 

Total temperature uncertainty was estimated at ± 0.25°C, due to sensor 
imprecision, ice point reference junction uncertainty, and temperature dis- 
tortions caused by the sensor in the measuring medium. 

Voltage and Current 

Uncertainties in the voltage and current were assumed to be .8 mV and 4mA 
based on statistical analysis of experimental results. 

Lengths 

The length and diameter of the test cylinder were assumed to have uncertain- 
ties of 5 and .03 mm, respectively. 

Pressures 

The uncertainty in pressure measurements were estimated based on manufac- 
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turer's accuracy statements and analysis of experimental results. 

Percent of Full Scale Uncertainty 

0 - 10% 1.2% of reading + 0.02% of full scale 

10 - 100% 1.4% of reading + 0.02% of full scale 

It was found that the maximum uncertainty in the Reynolds numbers and 

Nusselt numbers were 4.3% and 9.6%, respectively. The uncertainties at the 
different nominal Reynolds numbers are presented in Table 5. 


Table 5. Estimated uncertainties in the measured 
and Nusselt number. 

values of Reynolds number 

Reynolds number 

% Error in 
Reynolds number 

% Error in 
Nusselt number 

14,000 

2.8 

9.6 

9,000 

2.8 

8.4 

6,000 

2.8 

7.8 

4,000 

3.0 

7.4 

2,700 

3.4 

7.2 

1,800 

4.3 

7.0 


4. PRESENTATION AND DISCUSSION OF RESULTS 
In order to validate the experimental procedure, the convective heat 
transfer to the bare cylinder was measured in the wind tunnel facility over 
the Reynolds number range of 1000 to 20,000. The data were obtained using 

the specified constant heating rate of 1040 w/m 2 based on the outside sur- 
face of the test cylinder. Two sets of data were obtained in two "iden- 
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tical" wind tunnel facilities.* The measured Nusselt numbers as a function 
of Reynolds number are shown in Figure 12 along with the correlation curve 
developed by Morgan [7]. The data have a standard deviation of 6.1% 
relative to Morgan's correlation, thus validating the experimental pro- 
cedure. 

The maximum observed surface temperature variation (among the four 
thermocouples) for the bare cylinder tests was found to be 0.15°C. Later, 
it was found that for the shrouded cylinder tests, the maximum observed 
surface temperature variation was 0.3°C. These results validated the 
assumption that the surface temperature variation was not significant, thus, 
justifying utilization of only four thermocouples for surface temperature 
measurement. 

The convective heat transfer to shrouded cylinders was measured over 
the Reynolds number range of 1000 to 20,000 for three different sets of 
shrouds with radius ratios of 1.1, 1.4, and 2.1. Each set consisted of 
three shrouds with nominal ventilation factors of 0.09, 0.18, and 0.27, 
while the set of shrouds with a radius ratio of 1.1 had an additional shroud 
with a ventilation factor of 0.36. Each shroud was tested with two shroud 

orientations with respect to flow (e = 0°, e = 30°) as specified in Table 4. 

The set of shrouds with a radius ratio of 1.1 were also tested with shroud 

orientation of e = 15°. The results showing the variation of measured 
Nusselt numbers with Reynolds number and shroud orientation for each of the 
shrouds are shown in Figures 13 through 16. The correlation curve of Morgan 
[7] for convective heat transfer to bare cylinders is also shown in each 

figure. It should be noted that the Reynolds numbers were based on the bare 

cylinders diameter. 

*The wind tunnel was disassembled and moved to another building during the 
course of this investigation. 
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Fig. 13. (Continued). 
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Heat transfer data for perforated shrouds with ventilation factor of 0.36 (<j> = 1.1) 



The data indicate clearly that it is possible to enhance heat transfer 
to a solid cylinder by employing perforated shrouds (over the tested 
Reynolds number range). Convective heat transfer to the cylinder was en- 
hanced as high as 50% for the shrouded case compared to bare cylinder in 
some cases. It can be seen that shrouds with ventilation factor of .09 
retarded the convective heat transfer to the cylinder, while shrouds with 
ventilation factors of .18 and higher clearly enhanced the heat exchange. 
Shroud orientation influences the heat transfer with maximum heat transfer 
generally occurring when the holes are centered on either side of the stag- 
nation line (e = 30°), but that effect is of the second order compared to 
the ventilation factor and radius ratio. 

In order to interpret the ventilation factor more clearly, curves were 
fared through the experimental data for the shrouds with a radius ratio of 
1.1 for each shroud orientation over the Reynolds number range tested, 
recognizing the fact that the shroud ventilation factor has no effect at 
Reynolds numbers approaching zero. Then, the fared curves were used to 
estimate a "smoothed" Nusselt number for each ventilation factor and shroud 
orientation at selected Reynolds numbers between 1800 and 14000. Then the 
estimated Nusselt numbers were plotted against ventilation factor for dif- 
ferent Reynolds numbers, as shown in Figure 17a. In order to generate the 
curves shown in this figure two assumptions were made. First, it was 

assumed that Nusselt number should equal the bare cylinder Nusselt number at 
a ventilation factor of unity. Second, it was assumed that for a venti- 
lation factor of zero, the heat transfer is purely conductive between the 
shourd and the cylinder thus yielding a Nusselt number near unity. The 
Nusselt number represents the combined effects of heat conduction through 
the annulus and solid cylinder convection from the outside of the shroud. 
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Ventilation factor ( a ) 

Variation of extrapolated Nusselt number with ventilation factor for different Reynolds 
numbers, a) <|> = 1.1, b) $ = 1.4, and c) $ = 2.1. 
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Fig. 17. (Concluded). 




The shaded areas in Figure 17a represent the influence of shroud orienta- 
tion. From this figure it can be seen that for ventilation factors near 
0.2, the influence of the shroud on the convective heat transfer to the 
cylinder is favorable. Somewhat higher heat transfer occurs with the shroud 
in place than for the bare sensor at higher Reynolds numbers, but enhance- 
ment decreases as the Reynolds number decreases. For ventilation factors 
greater than approximately 0.2 the shroud becomes less effective. Further- 
more, for ventilation factors less than 0.1, the shroud retards convective 
heat transfer to the cylinder. The same type of behavior was generally 
observed for the shrouds with radius ratios of 1.4 and 2.1 as shown in 
Figures 17b and 17c. 

In order to study the combined effect of ventilation factor and radius 
ratio, it was decided to calculate an effective diameter, defined as the 
diameter which would yield a Nusselt number based on Morgan's correlation 
[7] equal to the experimentally measured Nusselt number, for each of the 
shrouded test cases. The results are presented in Table 6 where the average 
effective diameter of each shroud for the different shroud orientations and 
Reynolds numbers are presented as a function of ventilation factor and 
radius ratio. It was found that the effective diameter is basically inde- 
pendent of Reynolds number and shroud orientation. A dimensionless effec- 
tive diameter defined as the ratio of the effective diameter to bare cylin- 
der diameter is also presented in Table 6. The standard deviation of 
Nusselt numbers calculated using the given effective diameters in the 
Morgan's correlation with respect to the experimentally measured Nusselt 
numbers are also provided. The data indicate that the effective diameter 
could yield accurate estimations of the convective heat transfer for the 
perforated shroud case (in the range of experimental variables used). This 
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Table 6. Compilation of effective diameters for shrouds. 


Radius Ratio Ventilation 
Factor 


Effective 
Di ameter 
(mm) 


Standard deviation of Nusselt 
numbers calculated using the 
effective diameter in Morgan's 

Oimensionless correlation [7J with respect 

Effective to experimental data. (Values 

Diameter are in percentages). 


1.1 

.09 

5.68 ± .65 

.72 

5.65 

1.1 

.18 

9.41 ± 1.38 

1.20 

7.63 

1.1 

.27 

10.01 ± 1.31 

1.27 

7.08 

1.1 

.36 

11.08 ± .94 

1.40 

4.37 

1.4 

.09 

6.22 ± .46 

.79 

3.96 

1.4 

.18 

11.88 ± 1.05 

1.50 

4.86 

1.4 

.27 

14.75 ± 1.14 

1.87 

4.38 

2.1 

.09 

5.51 ± .54 

.70 

5.62 

2.1 

.18 

10.37 ± .78 

1.31 

3.95 

2.1 

.27 

12.47 ± .90 

1.58 

3.82 
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implied that the convective heat transfer for each of the shrouded cases 
differs from the convective heat transfer from the bare cylinder by the 
ratio of an effective diameter to the bare cylinders diameter raised to the 
appropriate power specified in Morgan's correlations [7]. Therefore, if the 
Reynolds nunber in Morgan's correlations is based on the given effective 
diameter, the thermal behavior of the shrouded cylinders can be predicted 
accurately (for the range of parameters specified here). 

The variation of the effective diameter with ventilation factor for 
different radius ratios is shown in Figure 18. Third order polynomial 
curves were fared through the experimental data for each radius ratio 
recognizing the fact that the dimensionless effective diameter should equal 
zero for a ventilation factor of zero, and unity for a ventilation factor of 
unity. It can be seen from Figure 18 that the relationship between 
effective diameter and ventilation factor for different radius ratios is 
consistent. The three third order polynomials were similar in shape and 
only varied from each other by constant factors. The data indicate that for 
ventilation factors between 0.2 and 0.4, the shrouds yield high effective 
diameters, thus enhancing convective heat transfer to the bare cylinder. 

For ventilation factors greater than 0.4, the shroud becomes less effective, 
while for ventilation factors less than 0.1, the shroud retards the convec- 
tive heat transfer. It can also be seen that the greatest enhancement of 
heat transfer occurs with a radius ratio of 1.4. Increasing the radius 
ratio above this value makes the shroud less effective. Another observation 
is that for a radius ratio of 1.4 with ventilation factors ranging from 0.2 

to 0.4, the convective heat transfer to the shrouded cylinder is higher than 
the convective heat transfer to a solid cylinder with a diameter equal to 
that of the shroud. 
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Variation of effective diameter with ventilation factor for different radius ratios 


The variation of effective diameter with radius ratio for different 
ventilation factors is shown in Figure 19. Third order polynomial curves 
were fared through the experimental data for each radius ratio using the 
requirement that the dimensionless effective diameter should equal unity as 
radius ratio approaches infinity. This requirement was based on the assump- 
tion that as radius ratio tends to infinity the effect of shroud on the 
convective heat transfer becomes negligible, given an adequate ventilation 
factor. It can be seen that the convective heat transfer to the cylinder is 
enhanced by using a radius ratio of around 1.4. 

Development of correlations between the effective diameter, ventilation 
factor, and radius ratio seemed inappropriate due to a limited number of 
data available. But the approximate variation of effective diameter with 
ventilation factor and radius ratio was obtained by incorporating the two- 
dimensional plots of Figures 18 and 19 into a three-dimensional plot as 
shown in Figure 20. It can be seen that the shroud enhances heat transfer 
to the bare cylinder for ventilation factors in the range 0.1 to 0.4 and 
radius ratio in the range 1.1 to 2.1. The maximum enhancement occurs for a 
radius ratio of 1.4 and ventilation factor between 0.2 and 0.3. 

5. CONCLUSION 

The influence of perforated cylindrical shrouds on the convective heat 
transfer to circular cylinders in transverse flow was measured experimental- 
ly. It was found that it is possible to estimate the convective heat trans- 
fer to a shrouded cylinder. The experiments showed that perforated shrouds 
with ventilation factors in the range 0.1 to 0.4 and radius ratios between 
1.1 to 2.1 could enhance the convective heat transfer to bare cylinders 
considerably. This enhancement in some cases produced heat transfer coeffi- 
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Fig. 19. Variation of effective diameter with radius ratio for different ventilation factors. 
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Fig. 20. Variation of effective diameter with ventilation factor and radius ratio. 



cients that were higher than those achieved by a cylinder having a diameter 
equal to that of the shroud. The maximum enhancement occurred with a radius 
of 1.4 and ventilation factor in the range of 0.2 to 0.3. 

The influence of radius ratio and ventilation factor appears to be 
systematic and predictable. It was found that the convective heat transfer 
for each of the shrouded cases differed from the convective heat transfer 
for the bare cylinder by the ratio of an effective diameter to the bare 
cylinder diameter raised to the appropriate power specified in Morgan's 
correlations [7], It was found that the effective diameters were indepen- 
dent of Reynolds number and shroud orientation and depended primarily on 
ventilation factor. 

Based on the limiting cases for the ventilation factor (0 and 1) and 
radius ratio (1 and infinity), it can be concluded that the experiments 
performed in this investigation were taken in the region where effective 
diameter changes the most. That is to say that the edges of the graph shown 
in Figure 20 can be drawn purely on the basis of physical logic. The dashed 
lines have been extrapolated based on that logic too. The fact that the 
surface described by the solid lines represents the dome of maximum influ- 
ence, suggests that effective diameter estimations can be made in the dashed 
region. 

Shroud orientation influences the heat transfer, but that affect is 
modest compared to the affects of ventilation factor and radius ratio. 

Effect of shroud orientation was more significant at higher Reynolds numbers 
and for ventilation factors between 0.1 and 0.3. The maximum heat transfer 
generally occurred when the shroud holes were centered on either side of the 
stagnation line. 
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APPENDIX A 


LONGITUDINAL TEMPERATURE VARIATION OF TEST CYLINDER 
One of the reasons for the wide dispersion in the published experi- 
mental results for convective heat transfer to bare cylinders can be attri- 
buted to heat conduction to the supports [7]. For negligible conduction 
losses, the aspect ratio (L/D) should be sufficiently large. Morgan [7] has 
reported that for a cylinder without heated end guard sections, the temper- 
ature is uniform over at least the center third of the cylinder when L/D > 
200, but that in some experimental investigations this ratio has been less 
than 10. Since the aspect ratio of the test cylinder used in this experi- 
ment was 26, the validity of assuming uniform temperature over the center 
third of the cylinder had to be verified. 

Knowing that the circumferential temperature variation over the test 
cylinder was negligible (Appendix B) , a two-dimensional finite volume nuner- 
ical approach [22] was used to calculate the steady state temperature dis- 
tribution of the test cylinder. Because of synmetry only half of the length 
of the cylinder was analyzed. The stainless steel support was assumed to be 
a cylinder with the same diameter as the test cylinder, exchanging heat with 
the surroundings along its circumference, and having its temperature at one 
end maintained at the free stream temperature. The geometry was discretized 
into an array of 8 radial and 13 axial nodes, as shown in Figure Al. Radial 
and axial conductances were calculated [22] using: 

C = kA/t (A.l) 

where 

k = thermal conductivity 
A = area of conduction 
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L = Length of test cylinder 


Fig. Al. Discretization of test cylinder for longitudinal temperature distribution analysis. 

cn 


l = length between nodes 


The interface conductivities were calculated using 


k = 


2k i k i+l 


k i + k i + l 


(A. 2) 


since the interfaces were placed midway between nodes [22], 

The thermal problem was solved using the Gauss-Seidel iteration with 
successive over relaxation [23] for both the lowest and highest Reynolds 
numbers used in the experiment -- 1,800 and 14,000. The longitudinal vari- 
ation of nondimensionalized surface temperature with respect to the temper- 
ature gradient between the average surface temperature and free stream tem- 
perature for the two cases under consideration are shown in Figure A2. It 
can be seen that the surface temperature along the center third of the 
cylinder is basically uniform. A listing of the utlized program is included 
here. 
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v>! 

ro 


PROGRAM 12 74/17 0 PT= 1* P OUND* A/ S/ M/-D*-DS FTN 5.1+587 * 86/0^/11. 13.47.44 

DU=-LJNG/-Ul»ARG a CUMM0N/-F1XED*CS“ USER/-FIXED*CB»-TB/-SB/-SL/-EP/-ID/-PMD/-ST*PL-50000 
FTM 5 » 1 "T31 # PL , 


1 

2 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
lb 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


PROGRAM T2 (INPUT* OUT PUT* TAPE 3 -OUT PUT ) 

DIMENSION T(0:i3*0:9)*XYl(12#B),XY2(12»8),ERRDP(12*e),TCm, 
1CR(12*0«6)»CZ(&*0«12)*C(12*8)*TCC(8)»R(0*9)»Z(13)»DZ(12)»ATC(8) 
M = 12 
N»8 
C 
C 
L 
C 

C THIS PROGRAM WAS USED FOR CALCULATING THE STEADY STATE 
C TEMPERATURE DISTRIBUTION UF THE TEST CYLINDER* NEGLECTING 
C CIRCUMFERENTIAL TEMPERATURE VARIATION 
C THIS PROGRAM CONSISTS OF THE FOLLOWING STEPS * 

C -SPECIFICATION OF RADIAL L AXIAL NODE LOCATIONS 
C -CALCULATION OF INTERFACE THERMAL CONDUCTIVITIES 
C -CALCULATION OF RADIAL £ AXIAL CONDUCTANCES 
C -SOLUTION OF PROBLEM USING GAUSS-SEIDEL ITERATION WITH 
C SUCCESIVE OVER-RELAXATION 
C 

C PARAMETERS USED HERE ARE: 

C RW = RADIUS OF HEATER WIRE 
C kill - RADIAL COORDINATES OF NODES 
C Z(J) * AXIAL COORDINATES OF NODES 
C 0 Z < J ) = AXIAL NODAL DISTANCING 
C W « RELAXATION PRAMETER 
C El = CONVERGENCE CRITERIA 
C OL = HEAT GENERATION PER UNIT LENGTH 


28 

29 

30 

31 

32 

33 

34 
3a 
3o 

37 

38 

39 

40 


C CK1 - THERMAL CONDUCTIVITY OF HEATER WIPE 

C CK2 « THERMAL CONDUCTIVITY OF POTTING COMPOUND 

C CK3 = THERMAL CONDUCTIVITY OF BRASS TUBE 

C CK4 ■= THERMAL CONDUCTIVITY OF STAINLESS STEEL SUPPORT 

C CH =■ C0VECT1VE HEAT TRANSFER COEFFICIENT 
C 

C TC(I) = RADIAL THERMAL CONDUCTIVITIES 

C TCC(J) - AXIAL THERMAL CONDUCTIVITIES 

C ATC(J) = AXIAL THERMAL CONDUCTIVITIES AT TUBE-SUPPORT 

C I NT EPF ACE 

C C R ( I * J ) = RADIAL CONDUCTANCES 

C CZ(I,J) - AXIAL CONDUCTANCES 

C C(I,J) » NODAL HLAT GENERATION TERMS 



PkOGkAM 1? 


74/17 


OPT = 1 » ROUND® A/ S/ M/-D,-DS 


FIN 5.1+587 


86/04/11. 13.47.44 


41 

42 

43 

44 
4 5 
4b 

47 

48 

49 

50 

51 

52 
33 

54 

55 
5b 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 
7o 

77 

78 

79 

--I 80 

0J 81 
82 


C f(l,J) - NODAL TEMPERATURE (DIFFERENCE OF NODAL £ FREE 
t STREAM r EMPERATURES ) 

C 

C 

C 

C SPECIFICATION OF SOME PARAMETERS 
C 

PI-4. *ATAN( 1.0) 

CK 1 = 1 9 • 

CK2-0. 52 
CK3-111. 

CK4= 14 . 4 
W- 1 . 5 
El-l.OE-5 
QL-26.0 
C 

X-25.4E-3 
RW» .020*X 
U1-0.010+X 
D2-0.021+X 
D3-.404+X 
04- . 3 99* X 
C H - 7 0 . 

C 

c 

c SPECIFICATION OF RADIAL NODAL COORDINATES 
C 

R ( 0 ) = 0 . 

R C 1 ) -0.0 
R ( 2 ) - D 1 

R(3)®R(2)+2.0*D1 
DO 1 1-4,6 
R ( I )»R ( 1-1 ) +D2 
1 CONTINUE 

R(7)-R(6)+2.0*U2 
P. (8 >-R (7)+D2 
P(9)-P (8) 

C 

C 

C SPECIFICATION OF AXIAL NODAL COORDINATES 

C 


Z ( 1 ) = 0 



•*vl 

4=» 


PROGK Art T 2 


74/1? U PT = 1* ROUND* A / S/ M/-D*-DS 


FTN 5.1+567 


86/04/11. 13.47.44 


83 

84 
65 


DU 2 1-2*11 
2 ( I ) - 2 ( I — 1 ) +D3 
CONTINUE 


66 

87 

86 

o9 

90 

91 

92 

93 

94 

95 

96 

97 
9b 
99 

i.00 

101 

102 

103 

104 

105 
i 06 
107 
1 ub 
109 
llu 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 


DO 201 1-12*13 
Z(I )«Z(I-1)+D4 
201 CONTINUE 
C 
C 

C CALCULATION OF AXIAL NODAL SPACING 

C 

C 

UZtl)=0.5*(Z(2)-Z(l) ) 

DO 3 1=2*10 
DZ(I>-D3 

3 CONTINUE 

DZ ( 11 ) - ( D3 + D4 ) / 2 . 

DZ { 12 ) =04 
C 
C 

C CALCULATION OF AXIAL £ RADIAL THERMAL CONDUCTIVITIES 

C 

C 

AA1-2 »0*CK1*CK4/(CK1 + CK4) 
AA2=2.0*CK2*CK4/(CK2+CK4) 

AA3»2.0*CK3*CK4/( CK3+CK4) 

DO 501 1-1,2 
ATC (I ) = AA1 

501 CONTINUE 

DO 502 1-3*6 
ATC { I ) • AA2 

502 CONTINUE 

00 503 1=7*8 
ATC ( I ) - AA3 
5C3 CONTINUE 
C 
C 

TC<1I*CK1 

TC ( 2) -2.0+CK1+CK2/ (CK1+CK2) 

DO 4 1=3*5 
TC ( I ) =CK2 

4 CONTINUE 

TC (fc)=2.0*CK2*CK3/(CK2+CK3) 



PROGRAM 12 


74/17 OP T»l* ROUND* A / S/ M/-D,-DS 


FTN 5.1+587 


86/04/11. 13.47.44 


125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 
15b 

159 

160 
161 
162 
163 

, 164 
<_n 165 

166 


TC ( 7) =CK 3 
C 

TCC ( 1 ) a C K 1 
TCC(2)=CK1 
DO 5 I m 3> 6 
TCC (1 )*CK2 

5 CONTINUE 
DO 6 I = 7 » 8 
TCC(I)*CK3 

6 CONTINUE 
C 

L 

C CALCULATION OF RADIAL CONDUCTANCES 

C 

C 

DO 7 1-1,10 

CR<1,8)«2.+PI*CH*R{8)*DZ(I) 

7 CONTINUE 

DO 8 1-1,10 
DO 9 J-1,7 

CR (I, J )»PI*TC(J )*DZ(1)*IR(J + 1)+R(J) ) / (R ( J +1 ) -R ( J I) 
9 CONTINUE 

8 CONTINUE 

DO 511 1-11,12 
OU 512 J * 1# 7 

CR ( I, J )=PI*CK4*DZ( I ) * IR (J+1)+R(J ) )/( R(J+1)-R CJ ) ) 

512 CONTINUE 
511 CONTINUE 

DO 513 I * 1 1 * 12 

CR( I, 8)=2.*PI*CH*R(8)*DZ( I) 

513 CONTINUE 
C 

C 

C 

C‘ CALCULATION OF AXIAL CONDUCTANCES 
C 


L 


DO 10 1-1,8 

B1«.25*PI*TCC (1 l/D3*( (P(II+R(I+1))**2.0-(RCH + P(I-1M**2.0) 
DO 11 J = 1, 9 
CZ(I»J )-Bl 
CONTINUE 
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P k 0 G k A H 12 

74/17- u PT *1 » R QUND* A / S/ M/-D.-DS FTN 5.1 + 587 86/04/11 

167 

10 

•CONTINUE 

16b 


DO 514 1*1,6 

169 


CZ ( I, 10)».25+P1*ATC < 1 1 / D34 ( (k< I)+R( I +1 )) **2 .0- C R ( I)+R( 1-111**2.01 

170 

514 

CONTINUE 

171 


00 51 5 1=1.6 

172 


DV=.25*PI*CK4/0 4*((R(I)+R(I + in**2.0-(Rm+RU-l)>**2.0) 

173 


DO 516 J=ll» 12 

174 


CZ ( It J ) =3V 

175 

bio 

CONTINUE 

176 

515 

CONTINUE 

177 

C 


176 

C 


179 

L 


180 

C CALCULATION OF NODAL HEAT GENERATION TERMS 

181 

C 


182 

C 


183 


DO 44 1=1. M 

184 


DO 45 J * 1 » h 

18b 


C(I.J)=0.0 

186 

45 

CONTINUE 

167 

44 

CONTINUE 

168 

C 


18V 


B3=0L*(0.5*R(2)/RW)**2.0 

190 


B4=QL*(1-(0.5*R(2)/RW)**2.0) 

191 


DO 13 1=1.10 

192 


C(I,1)=B3*DZ(I) 

193 


C ( 1.2 ) =U4*DZ ( I ) 

194 

l 3 

CONTINUE 

195 

C 


196 

C 


197 

C INITIALIZING THE PROBLEM FOR GAUSS-SEIDEL PROCEDURE 

198 

C 


199 


DO 14 1=1, M 

200 


CR ( I, 0 J«0. 

201 


T ( 1 , 0 ) ■ 1 . 

202 


T( 1,91=0. 

203 

1*. 

CONTINUE 

204 


DO 16 J-1,8 

205 


T ( 0, J )»1. 

206 


CZ( J,0»=0. 

207 


T ( 13, J ) =0. 

208 

15 

CONTINUE 



PROGRAM T 2 


74/17 


OPT *1»RQUND* A/ S/ M/-D.-DS 


FTN 5.1+507 


86/04/11. 13 .A7.<*<. 


c o9 

C 


21u 


00 10 I ■ 1 * M 

211 


DO 16 J» 1 » N 

212 


T(I»J)=0.0 

in 


XYU1»JJ-0.U 

ilk 

lb 

CONTINUE 

215 

lb 

CONTINUE 

216 

c 


217 

c 


2 IB 

t 


214 

c 

iOLUTICN OF PROBLEM USING GAUSS-SEIDEL ITERATION 

220 

c 

WITH SUCCESSIVE OVER-RELAXATION 

221 

c 


222 

l 


223 


DO 200 K= 1 » 50000 

224 

c 


226 


DO 20 I-l.M 

226 


DO 21 J-l.N 

22 7 


TCI»J)-IC(i#J)+CRII»J-l)*T (It J-l)«CR(If J)*T(I#J+11+CZ(J,I»* 

228 


iTII + l»J)+CZ(J>I-i)*T(I-l.jn/(CR(I»J-l)+CRII#J)+CZ(J#n + C2UfI 

229 


XY211# J)-T( 1 > J ) 

230 

21 

CUNTINUF 

2 31 

20 

CONTINUE 

232 

C 


233 


DO 25 I ' 1 » H 

234 


DO 26 J-l.N 

235 


ERR UR ( I»J)*A8S(XY2(I>J)-XY1(I»J) ) 

236 


IF ( ERROR ( I» J) .GT. El ) THEN 

237 


GO TO 71 

23b 


ELSE 

234 


IF ( I .EO.M.AND.J .EO.N) THEN 

240 


GO TO 70 

241 


ELSE 

242 


GO TO 26 

243 


END IF 

244 


END IF 

245 

26 

CONTINUE 

24b 

25 

CONTINUE 

c 47 

71 

CONTINUE 

24B 


DU 27 1*1. N 

249 


DO 28 J* 1. N 

250 


XYHI»J)*XYl(l.J)+W+(XY2( I.JI-XY1II.J) ) 





00 


PRUGk AM T £ 


76/17 UPT-1, ROUND* A/ S/ M/-D,-DS F TN 5.1 + 567 


86/06/11. 13.67.66 


£51 


T ( I » J ) -XY1 { 1, J ) 

Z-3l 

£8 

CGNTINUr 

l 53 

£7 

CONTINUE 

£5*i 

2C0 

CONTINUE 

£55 


HR ITE ( 3. 601 

£56 

60 

FORMAT ( 3 X » ' MOT CONVERGING*) 

£57 


GO TO 00 

£ ob 

70 

CUNTINUT 

£ 59 


DU 29 1*1, M 

£ 60 


DO 30 J- 1, N 

£61 


WRITE(3,6l) 1 , J , T ( I » J ) 

2 62 

61 

FORMAT (3X, «T(»» 12#', »#12> • »-■ 

£63 

30 

CONTI NUE 

£ 66 

2 9 

CONTINUE 

£65 

80 

STOP 

£ 6b 


END 



APPENDIX B 


CIRCUMFERENTIAL TEMPERATURE VARIATION OF TEST CYLINDER 

Local heat transfer coefficient on a circular cylinder varies consider- 
ably across the surface of the cylinder [9]. This would in turn cause cir- 
cumferential variation of surface temperature around the cylinder. Even 
though a relatively high thermal conductivity tube (brass) was utilized to 
prevent significant surface temperature variation, further analytical inves- 
tigation was required to determine maximum steady state surface temperature 
variation for the bare cylinder. A two-dimensional finite volume numerical 
analysis [22] of the test cylinder was utilized for this purpose. The geom- 
etry was discretized into an array of 8 radial nodes and 7 circumferential 
nodes, as shown in Figure Bl. The radial and circumferential conductances, 
and the interface thermal conductivities were calculated as prescribed in 
Appendix A. Schmidt and Wenner's [9] data for the circumferential variation 
of Nusselt number for a Reynolds number of 39,800, shown in Figure 10, was 
used in this study. 

The thermal problem was solved using the Gauss-Seidel iteration with 
successive over relaxation [23]. It was found that the maximum circumfer- 
ential temperature variation would be 0.03°C for the considered case. A 
listing of the utilized program is given here. 
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Fig. Bl. Discretization of test cylinder for circumferential temperature distribution analysis. 


PROGRAM li 7A / 17 OPT *1,PUUND« A / S/ M/-D»-DS UN 0.1 + 587 86/05/01 . 19.52.3A 

uj*-U!NG/-u I mFG* Ct,MNnN/-E IXtD.CS* US L R / -E 1 X 1 0, OB «-l B/ -S B / -S L /-E R/- I D /-PMD /~S T * PL » 50000 
* IN3*1*I1*Pl. 


1 

2 

3 

5 

■j 

o 

7 

d 

•y 

id 

11 

12 
Id 
iH 

1 3 

io 

17 

iu 
1 + 
20 
21 
22 
4 d 

2 A 
2 3 
26 
27 
2 o 
20 
dO 
dl 
d2 
d j 
3a 
33 
d 6 

d 7 
3d 


Pk LIon A.1 T1 ( INruI » uU T PUT * TAPE 3- GUT PUT ) 

D I Ml NS I I'll T(d,7)»XYl(b»7),>Y2(8»7)»£RR0R(8*7)*TC(6>»TCC<8)»R<8) 
1 ,KC (<J,7)»CC 18* 7)#i.(b# 7) #CH(7» 

11*8 
N= 7 

L 

L 

C 

L 

C IhlS PROGRAM WAS UStl) EUR CALCULATING THE STEADY STATE 

C It MPEKATDRL distribution of the test cylinder* neglecting 

C AXIAL TENPERATUKt VARIATIUN 

L I H I S PROGRAM CONSISTS OF THE FOLLOWING STEPS: 

C -SPECIFICATION UF RADIAL L CIRCUMFERENTIAL NODE LOCATIONS 
L -CALCULATION OF INTtREACfc THERMAL CONDUCTIVITIES 
L -CALCULATION OF RACIAL l CIRCUMFERENTIAL CONDUCTANCES 
C -SOLUTIUN UP PRUBLtfl USING GAUSS-SEIOEL ITERATION WITH 
C SUCCESUF OVER-RlL AXATIUN 
C 

C PAPAMtTtRS OSH HE R t ARE: 

L h» ■ RADIUS OF HtATtK WIRE 
C R ( 1 ) = RADIAL COORDINATES OF NODES 
C Cl * C IRCUMT EPtNT 1AL DISTANCING 
C * * RELAXATION PARAMETER 
Cii* CONVERGENCE CRITERIA 
C .L * HEAT GENERATION PER UNIT LENGTH 

C CRl * THERMAL CONDUCTIVITY UF HEATER WIRE 

C Cn2 = THERMAL CONDUCTIVITY QE POTTING COMPOUND 

c C R 3 * THERMAL CONDUCTIVITY UE BRASS TUBE 

C Kill * R Alt 1 AL THERMAL CONDUCTIVITIES 
C TCCII) = CIRCUMFtRLNTIAL THERMAL CONDUCTIVITIES 
L iHUI » CONVECTIVE HEAT TRANSFER COEFFICIENTS 
C R C ( 1 « J ) = RADIAL CONDUCTANCES 
C LC(I.J) = C IRCOMFEPENTIAL CONDUCTANCES 
C C(I.J) = NODAL HEAT GENERATION TERMS 

L l(l,j) = NODAL ItMFtRATUFES (DIFFERENCE OF NODAL E FREE 
C STREAM TEMPERATURES) 

C 


j v C 

AG L 


CO 



CO 

ro 

PkOopAM fi /h/17 LPI *1# ROUND* a/ S/ M/-L»> -DS FTN 5.14587 86/05/01 . 19.52.3S 


*1 i 

t j 

SS 
S5 
So 
S 7 
sa 

S'V 
> J 

l' A 

32 
3 J 


L 

C 

L . i * CIFICAT 1IJN tU HARAHlTFFS 
C 
( 

FI = 4.4A|AM1.C» 

DWI/c.O 
( M * 1 9 . 

C K 2 = 0 . 3 2 
CKj=m. 
c - 1 . j 
E 1 = 1 . 01 -b 
01 =26.0 


34 
33 
3b 
37 
3 o 

5V 

UO 

cl 
u2 
O 3 
o'* 
o5 
tfa 
c 7 
o J 
67 
7c 
71 
7 c 

73 

74 
7 3 
7c 
77 
7 d 
77 
BC 
fal 
02 


0 

7=25.41 -3 
R*=.U2C*X 
Di=O.OiO*X 
D2 = u.c<:1*7. 

L 

C ji-ECI F1C AT UlN UP KA01AL NODAL CUORDlNAT ES 
0 

R 111 = G . 0 
ft < 2 ) = l'l 

kl3J«ftC2l+2 .0*01 
DU IOC I = 4# 6 
ft< I)*ft. I 1-1) +02 
UH CONTINUE 

M7)*Plt-)*2.0*D2 
ft m = K(?)+o? 
c 
c 

C SrtCiPiCAT I UM Of CONVECTIVE HEAT TRANSFER COEFFICIENTS 
c FROM SCHMIDT l RlNNEH'S LATA 191 

C 

C 1 1 < 1 ) = 6 3 3 . 

CHI2) *540. 

CHIB) *343. 

C H ( 4 > = 164. 

CM I 3 I * J 2 7. 

C F { o I = 4 5 b . 

CM 7) =055. 

C 



h K O k> K A 11 


74/17 


uHT>l»RUUm)« A/ S/ M/-D,-DS 


FTN 5. 1*587 


86/05/01. 19.52.34 


cl 
d** 
b j 
a i> 
o 7 
bo 
b9 

•/O 

91 

92 

93 

94 
93 
9o 

97 

98 

99 
. Ob 
l 01 

* u2 

1 04 
1 03 
1 UO 

x o 7 

I JH 
10 9 

I I J 

III 
lit 

1 13 
i 14 

1*3 

lie 
117 
1 lb 
i 19 
1 1 2u 

00 i ltl 

w i 122 

i 1 2d 

* 24 


C 

0 oalculati oti ;ir i-auial t h t mi a l conductivities 

c 

t 1 1 2 1 = a i 

TL(3l=2.0»LKi*oKt/lCKl+CK2) 

UL! 1 I * 4 , o 

TC(1I*CK2 

1 CUNTINUI 

TO (7)=2.0*LK2»CK3/ (CK2+CK3) 
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APPENDIX C 


RADIATION CORRECTION FOR BARE CYLINDER 
The bare cylinder exchanges radiant energy with the wind tunnel walls 
as shown in Figure C.l. The bare cylinder is designated as surface 1, and 
wind tunnel walls are designated as surfaces 2 and 4. Surfaces 3 and 4 are 
imaginary surfaces -- at free-stream temperature -- for completing the 
enclosure. The formulation is based on the assumption of radiant exchange 
between gray surfaces through non-participating media. 

The net radiant heat flux for each of the surfaces can be obtained 

from: 


<q,) n et = >>, - Si (C.l) 

The irradiation, G, on each surface is due to energy leaving that surface 
(if it can "see" itself) and the other surfaces: 


A.G. = 7 A. F. . 
11 j J J - n 


Applying reciprocity to this equation yields: 


(C.2) 


G. = l F. . J. 
1 j ^ J 


(C.3) 


The radiosity, J, is the sum of the emitted energy and reflected irradi 
ation, and can be calculated according to: 


J. 

i 


£ . 
1 



(1 - e.) 


G. 


(C.4) 
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Aluminum wall with enamel paint 



Aluminum wall 

Fig. Cl. Schematic diagram for radiation heat flux correction 
for bare cylinder case. ^erection 
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Combining Equations C.l, C.3, and C.4 results in 


(d^net ” " ( J i ” E b,i^ ^.5) 

1 - e i 

The radiosity for surface 1 can be calculated according to 

Jl = e ct + ( 1 -e) [Fj-2 J2 + F1-3 1 > R3 ' + F1-4 R4] (C.6) 

where surface 3' is the combination of surfaces 3 and 5. 

The radiosities for the other surfaces are: 

J 3 ' = qTj* (black body assumption) (C.7) 

J 2 = + ( l-e 2 ) [F 2 -iJi + F 2 _ 3 ’ J 3 ' + F 2 _4 J 4 ] (C. 8 ) 

J 4 = £ 4 ^ 1 ^ + (1-e^) [ F 4 -1 + F 4 _2 J 2 + F 4 _3 1 J 3 1 ] (C.9) 

Combining Equations C.5 through C.9 results in 

q Rl = ^Pnet = " 0>562 CT f°* 430 + °- 215 T w " °* 639 T4 S ] (C.10) 

where, crossed-string method [ 12 ] has been used to calculate the shape 
factors (F. .), and the following total normal emissivities have been used: 

* ”J 


e = .36 

e 2 = 0.9 (enamel paint) [12] 

£3 = 0.07 (aluminum, rough plate) [12] 

Equation C.10 was used for calculating the radiation heat flux corrections 
for the bare cylinder case. 
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